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Abstract The Helioseismic and Magnetic Imager (HMI) instrument onboard the Solar Dy- 
namics Observatory (SDO) satellite is designed to produce high-resolution Doppler- velocity 
maps of oscillations at the solar surface with high temporal cadence. To take advantage of 
these high-quality oscillation data, a time - distance helioseismology pipeline (Zhao et al . , 
Solar Phys. submitted, 2010) has been implemented at the Joint Science Operations Center 
(JSOC) at Stanford University. The aim of this pipeline is to generate maps of acoustic travel 
times from oscillations on the solar surface, and to infer subsurface 3D flow velocities and 
sound-speed perturbations. The wave travel times are measured from cross-covariances of 
the observed solar oscillation signals. For implementation into the pipeline we have inves- 
tigated three different travel-time definitions developed in time - distance helioseismology: 
a Gabor-wavelet fitting (Kosovichev and Duvall, SCORE’96: Solar Convection and Oscil- 
lations and Their Relationship, ASSL, Dordrecht, 241, 1997), a minimization relative to a 
reference cross-covariance function (Gizon and Birch, Astrophys. J. 571, 966, 2002), and 
a linearized version of the minimization method (Gizon and Birch, Astrophys. J. 614, 472, 
2004). Using Doppler-velocity data from the Michelson Doppler Imager (MDI) instrument 
onboard SOHO, we tested and compared these definitions for the mean and difference travel- 
time perturbations measured from reciprocal signals. Although all three procedures return 
similar travel times in a quiet-Sun region, the method of Gizon and Birch ( Astrophys . J. 614, 
472, 2004) gives travel times that are significantly different from the others in a magnetic 
(active) region. Thus, for the pipeline implementation we chose the procedures of Koso- 
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vichev and Duvall ( SCORE’ 96 : Solar Convection and Oscillations and Their Relationship, 
ASSL, Dordrecht, 241, 1997) and Gizon and Birch ( Astrophys . J. 571, 966, 2002). We in- 
vestigated the relationships among these three travel-time definitions, their sensitivities to 
fitting parameters, and estimated the random errors that they produce. 

Keywords Sun: helioseismology • Sun: time - distance analysis • HMI 


1. Introduction 

Time - distance helioseismology (Duvall et al, 1993) allows access to the subsurface 
physical properties of the quiet Sun (see e.g. Kosovichev and Duvall, 1997), includ- 
ing supergranulation (see e.g. Duvall and Gizon, 2000; Zhao and Kosovichev, 2003; 
Hirzberger et al, 2007), and also the subsurface properties of active regions (see e.g. Koso- 
vichev, Duvall, and Scherrer, 2000; Zhao, Kosovichev, and Duvall, 2001; Jensen etal, 2001; 
Couvidat, Birch, and Kosovichev, 2006). So far the main goal has been to determine the ve- 
locity of material flows in both the horizontal and vertical directions, and the sound-speed 
profile as a function of depth. As such, time - distance helioseismology is a major tool to 
achieve the primary science goal of the Helioseismic and Magnetic Imager (HMI) instru- 
ment onboard the Solar Dynamics Observatory (SDO) satellite: to study the origin of solar 
variability, and to characterize and understand the Sun’s interior and various components of 
magnetic activity. Therefore, a time - distance helioseismology data-analysis pipeline (de- 
scribed by Zhao et al. , 2010) has been implemented by the authors, at the SDO Joint Science 
Operations Center (JSOC). It will process the line-of-sight Doppler- velocity (and potentially 
continuum intensity) maps obtained by HMI to generate wave travel-time maps of the solar 
surface and invert these maps to generate estimates of the flow velocities and sound-speed 
perturbations beneath the surface. 

The travel times of acoustic or surface-gravity wavepackets are measured from tempo- 
ral cross-covariances of solar-oscillation signals measured at different positions on the solar 
disk. Over the years, various definitions of these travel times have been designed (Koso- 
vichev and Duvall, 1997 [hereafter KD97]; Gizon and Birch, 2002 [hereafter GB02]; Gi- 
zon and Birch, 2004 [hereafter GB04]). They generally return different times for a given 
cross-covariance function. Upon their implementation in the HMI time - distance pipeline 
we compare the results of these definitions using Doppler- velocity data from the Michelson 
Doppler Imager (MDI) instrument (Scherrer et al., 1995). We compare these three travel- 
time definitions for both quiet-Sun and active -region data. We only compute the mean and 
difference travel-time perturbations in the center-to-annulus geometric scheme, neglecting 
the North - South and East- West travel-time differences (see Couvidat and Birch, 2009, for 
a comparison of these travel-time differences in the presence of artificial uniform and steady 
flows, and for the three travel-time definitions). Moreover, we only consider acoustic waves 
(p modes). Jackiewicz et al. (2007) already studied the travel times of surface-gravity waves 
(/ modes). Finally, Roth, Gizon, and Beck (2007) compared the travel-time differences re- 
turned by the KD97 and GB04 methods in the quiet Sun. In this paper we expand part of 
their analysis to include the GB02 approach, to include mean travel-time perturbations, and 
we also work on active -region data. Other possible ways to measure the travel-time shifts 
of wavepackets are not considered here. For instance, it is possible to compute the instan- 
taneous phase of the cross-covariance using the Hilbert transform, and to define the travel 
time as the time of zero-crossing of this instantaneous phase (see e.g. Duvall et al, 1996). 
Other related techniques are phase-correlation holography (Findsey and Braun, 2005) and 
acoustic imaging (Chou and Duvall, 2000). 
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In Section 2 we review the time - distance formalism and discuss the travel-time defin- 
itions. In Section 3 we present the implementation of these travel-time measurement pro- 
cedures in the HMI JSOC data-analysis pipeline. In Section 4 we compare the results of 
these different approaches for quiet-Sun and active-region data from SOHO/MDI. Finally, 
we present our conclusions in Section 5. 


2. Time - Distance Formalism 

2.1. Computation of the Cross-Covariances of Solar Oscillations 

In this section we outline the cross-covariance calculations used for testing our travel-time 
measurement procedures and codes. We apply the measurement procedure detailed by Cou- 
vidat, Birch, and Kosovichev (2006). The measurements of the travel times are based on 
fitting the cross-covariances [C(y\,Y 2 , t) \ between source | r 1 1 and receiver [ r 2 1 calculated 
from a line-of-sight-velocity data cube [0(r, t) ] where r is the horizontal position vector, 
and T is time: 


Here T is the temporal duration of the observation. A high-pass filter is applied at 1.7 mHz 
to the Fourier transform of the data cube, to discard supergranulation. An /-mode filter is 
applied to discard the /-mode signal. To reduce the noise level of the cross-covariances, 
<p{ r, t) is filtered in the Fourier domain using a Gaussian phase-speed filter (Duvall et al., 
1997) for short travel distance A = |m — it|. Then, the point-to-point cross-covariances 
are averaged over annuli of radius A centered on the source. Such point-to-annulus cross- 
covariances are computed for 55 distances [A] ranging from 3.7 Mm to 66.7 Mm. Only 11 
phase-speed filters are applied for these 55 distances (one filter for five distances). These 
distances are not those selected for the time - distance pipeline: they are used here only to 
test our implementation of the travel-time definitions. To further increase the signal-to-noise 
ratio for the point-to-annulus cross-covariances, the latter are shifted in time - so that they 
are centered on specific reference times - and averaged in groups of five distances (the five 
A computed with the same phase-speed filter), thus giving only 11 final A values for mea- 
suring the acoustic travel times. For instance, five cross-covariances computed at distances 
ranging from 3.7 to 8.7 Mm are shifted in time to be centered on the same reference time 
(here the group travel time of a wavepacket propagating in the quiet Sun over the distance 
A = 6.2 Mm, as measured by the Gabor- wavelet fit presented in the next section) and are 
averaged to produce the final point-to-annulus cross-covariance at A = 6.2 Mm. The values 
of A and properties of the phase-speed filters used in this paper are summarized in Table 1. 

2.2. Travel-Time Definitions Implemented in the Pipeline 

The travel times of the acoustic wavepackets are measured from the point-to-annulus cross- 
covariances. Historically, the cross-covariances were fitted by a Gabor wavelet defined as 
(Kosovichev and Duvall, 1997) [KD97]: 



( 1 ) 
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Table 1 Parameters used for the time - distance analysis: central source - receiver distances [A], range of 
A, properties of the Gaussian phase-speed filters applied to the Fourier transform of the line-of-sight- velocity 
data cube (no is the central phase speed and FWHM is the full width at half maximum), and centers [to] of 
the temporal window applied to isolate the first-bounce skip on the cross-covariances. 


A (Mm) 

range (Mm) 

no (km s 1 ) 

FWHM (km s” 1 ) 

to (min) 

6.20 

3.70-8.70 

12.77 

6.18 

19.0 

8.70 

6.20-11.2 

14.87 

6.18 

23.3 

11.60 

8.70-14.5 

17.49 

6.18 

24.4 

16.95 

14.5-19.4 

24.82 

9.09 

28.7 

24.35 

19.4-29.3 

35.46 

12.36 

33.5 

30.55 

26.0-35.1 

39.71 

7.18 

36.3 

36.75 

31.8-41.7 

43.29 

7.42 

38.7 

42.95 

38.4-47.5 

47.67 

8.41 

40.8 

49.15 

44.2-54.1 

52.26 

10.5 

42.8 

55.35 

50.8-59.9 

57.16 

8.90 

44.7 

61.65 

56.6-66.7 

61.13 

8.03 

45.0 


where A is the amplitude of the wavelet, coq is the central temporal frequency of the 
wavepacket, Sco is a measure of its frequency width, r p is the phase-travel time, and r g is the 
group travel time. Commonly, only phase-travel times are used for the time - distance inver- 
sions, because they can usually be determined more accurately than the group travel times. 
However, r p is not unique and is defined modulo 2tt/coq. Therefore, the following rule has 
been used in this paper to select the phase-travel times: the value of r p that is selected is the 
closest to - and if two are equally distant, smaller than - the group travel time r g . 

Subsequently, two definitions of the wavepacket travel times based on studies in geo- 
physics were added. Gizon and Birch (2002) [GB02] define the travel time r(r, A) as: 

r±(r, A) = argmin X±(r, A, t) (3) 

t 


where A ± (r, A, t) are the functions: 

X±(r, A, t) = f(±t')[C( r, A, t’) - C ref (A, t’ T 0] 2 (4) 

t' 


and f(t) is a one-sided window function that selects the first-bounce skip of the cross- 
covariance (the skip corresponding to waves traveling directly from the center to the annu- 
lus) for the positive-time (i.e. outgoing waves) or negative-time (i.e. ingoing waves) branches 
of the cross-covariance. The same window function is used when fitting the cross-covariance 
by a Gabor wavelet even though it does not explicitly appear in Equation (2). The function 
C ref (A, t) is a reference cross-covariance. Here, C ref (A, t) is calculated as a spatial average 
of the measured cross-covariances over a region of the quiet Sun. 

Lastly, Gizon and Birch (2004) [GB04] proposed another definition of the travel time 
that is linear with respect to the cross-covariance: 

r±(r, A) = J2 W±(A, 0[C(r, A, t) — C ref (A, t)] (5) 
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where W± is a weight function: 


W ± (A,t) = 


Tf (±QC ref (A, t) 

j:,, f(±t')[c^(A,nr- 


(6) 


where C ref ( A , t) denotes the temporal derivative of the reference cross-covariance. Essen- 
tially, GB04 is a linearization of GB02 with respect to a parameter e when e — > 0, where e 
is used to define a smooth cross-covariance C f (r, A, t): 


C e (r, A, t) = eC(r, A, t) + (1 - e)C ref (A, t) 


(7) 


The definition GB02 applied to C e (r, A ,t) simplifies to GB04 in the limit e — > 0. For 
noiseless data, i.e. for an observation time T = +oc, and when the amplitude of the cross- 
covariance does not vary, GB02 and GB04 are equivalent. GB04 can be applied only when 
the amplitude variations of C(r, A , t) are small. This is not the case in active regions. While 
the travel-time variations in sunspots do not exceed, roughly, 20% of the typical wave pe- 
riod, the cross-covariance amplitude changes by more than 100%. With GB04 there is a 
linear relationship between the travel times and C(r, A, t) — C ref ( A, /). 

As emphasized by the use of the ± symbol, for the three definitions of the travel times 
we measure r separately for the outgoing (r + ) and ingoing (r_) waves. We define the 
mean travel time, r mean (r, A) = (r+(r, A) + r_(r,A))/2, and the difference travel time, 
Tdiff(r, A) = r + (r, A) — r_(r, A). The corresponding travel-time perturbations Sr are pro- 
duced by subtracting the travel time obtained from the reference cross-covariance. For in- 
stance, <5r mean (r, A) = r mean (r, A) - T mean;r ef(A) where t mean;ref (A) is the mean travel time 
measured for the reference cross-covariance. 


3. Implementation of the Travel-Time Definitions in the HMI Time - Distance 

Pipeline 

The three travel-time definitions were implemented in the Fortran codes of the HMI 
time - distance helioseismology pipeline, based on the Data Resource Management Sys- 
tem (DRMS) at Stanford University. They represent subroutines of the main time - distance 
program, which is also written in Fortran and encapsulated inside a C-wrapper to access the 
DRMS I/O functions. The input parameters of the travel-time subroutines include the cross- 
covariances, and their output are the travel-time maps in FITS format. The mean and dif- 
ference travel-time maps are calculated, as well as the North - South and East - West differ- 
ence travel-time maps in a center-to-quadrant geometry (used to measure horizontal flows). 
The Gabor- wavelet fitting procedure (KD97) is a least-squares fit using the Fevenberg- 
Marquardt algorithm. The GB02 and GB04 procedures use a Whittaker- Shannon interpo- 
lation formula to interpolate the temporal cross-covariances from a grid with a 45-second 
sampling rate (nominal HMI cadence) or a 60-second sampling rate (MDI cadence) onto a 
grid with a five-second sampling rate. The interpolation subroutine computes interpolated 
points only in a narrow time range around the temporal window chosen to select the first- 
bounce skip (see, e.g., the last column of Table 1 for some typical values). This speeds up 
the computations. 

For the GB02 definition, the code uses the interpolated cross-covariances to find the 
value of t = train at which X±(r, A , t) is minimum. Following Gizon and Birch (2002) we 
then fit X±(r, A, t) with a parabola around r 11lin (using five measurement points). The code 
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returns an estimate of the parameters cq, c i, and c 2 of the parabola y = cy + C| / + ci_ t 2 , 
and the travel-time measurement is given by the location of the minimum of this parabola 
T± = Ci/(2c 2 ). 

For the GB04 definition, the time derivative C ref (A, t) is numerically computed for the 
interpolated reference cross-covariance using a five-point equation. 

The three travel-time definitions return outliers (misfits in the case of KD97 and GB02). 
To correct for these outliers, each travel-time subroutine computes the mean [/x] and stan- 
dard deviation [er] of the ingoing and outgoing travel times. When | r± (r, A) — /z± | is greater 
than X times o± (where A is a parameter that can be set in the code; its value is typi- 
cally larger than three), the cross-covariance at location r is averaged over the eight closest 
neighboring locations. Then the travel times at r are re-computed using the averaged cross- 
covariances. 


4. Comparison of the Travel-Time Definitions 

To test our implementation of the travel-time definitions, we apply them to MDI Doppler- 
velocity data. This allows a comparison of the travel-time definitions in the quiet Sun and in 
an active region. 

4.1. Quiet-Sun Region 

We work on a line-of-sight Doppler- velocity data cube obtained from MDI high-resolution 
(Hi-Res) data taken in May 2001 for a quiet-Sun region (Gizon and Birch, 2004). By quiet 
Sun we mean a region devoid of any significant magnetic activity (e.g . , no sunspot, no plage). 
The observation duration is 512 minutes (or 8.5 hours). We choose the Hi-Res data because 
the spatial sampling rate (0.6" per pixel) is close to the HMI sampling rate. The data were 
remapped onto the heliographic coordinates using Postel’s projection, and the region was 
tracked at the Carrington rate. After rebinning, the spatial resolution of our data cube is 
Sx = Sy = 0.826 Mm, and the temporal resolution is St = I minute. The dimensions of this 
data cube are 256 x 256 x 512 where the first two dimensions are the horizontal ones, and 
the last dimension is the temporal one. We apply the time - distance formalism on this data 
cube. By default, the temporal windows used to select the first-bounce skip on the cross- 
covariances are 20 minutes wide. Figures 1 and 2 show the mean and difference travel-time 
perturbation maps obtained for two distances A = 6.2 and A = 30.55 Mm using the three 
travel-time definitions. All of the maps look very similar: they exhibit the same features at 
the same locations. However, the amplitude of these features varies. At A = 6.2 Mm, the 
amplitude of <5tdiff(r, A) in the supergranules (black spots on the maps) is larger in absolute 
value for GB04 than for GB02 and KD97. Table 2 lists the values of the Pearson linear 
correlation coefficients between the different travel-time maps obtained for the quiet Sun. 
These coefficients are high, especially between GB02 and GB04. 

Figure 3 shows histograms of the mean and difference travel-time perturbations for 
A = 6.2 and A = 30.55 Mm. The travel-time distributions for the three definitions look very 
similar (more or less close to a Gaussian distribution). The standard deviations are larger for 
GB02 and GB04 than for KD97 at A = 6.2 Mm, but they are smaller at A = 30.55 Mm. 
Therefore, while GB02 and GB04 seem to produce slightly “noisier” travel-time maps at 
short distances than KD97, this is the opposite for larger distances (what matters is actu- 
ally not the standard deviation - an assessment of the noise level - but the signal-to-noise 
ratio). This is consistent with the results of Roth, Gizon, and Beck (2007) for very large 
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Figure 1 Mean (upper panels) and difference (lower panels) travel-time perturbations for the three travel-time definitions (left panels: KD97; central panels: GB02; right panels: 
GB04) and A = 6.2 Mm. 
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GB04) and A = 30.55 Mm. The scale has been truncated to f— 35, 35] seconds to highlight the features in the maps. 
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Table 2 Linear Pearson correlation coefficients between the different travel-time maps in the quiet Sun, for 
the entire 256 x 256 maps. 

Map 

KD97/GB02 

KD97/GB04 

GB02/GB04 

^ r mean A = 6.2 Mm 

0.781 

0.748 

0.956 

<5 r mea n A = 30.55 Mm 

0.783 

0.726 

0.940 

cSr^iff A = 6.2 Mm 

0.876 

0.839 

0.945 

5 r^iff A = 30.55 Mm 

0.849 

0.797 

0.913 

A = 6.2 

Mm 

A = 6.2 Mm 




-30-20-10 0 10 20 30 

<5T mean (seconds) 




-30-20-10 0 10 20 30 

(5r diff (seconds) 

A = 30.55 Mm 



-30-20-10 0 10 20 30 

<5 T mean (seconds) 


-30-20-10 0 10 20 30 

c5 T dif f (seconds) 


Figure 3 Histograms of the mean (left column) and difference (right column) travel-time perturbations for 
the two distances A = 6.2 (upper panels) and A = 30.55 (lower panels) Mm. Black lines are for GB02, green 
lines are for GB04, and red lines are for KD97. 


distances (24 heliocentric degrees, or about 291 Mm). For short distance, here A = 6.2 Mm, 
the distributions of the mean travel-time perturbations are noticeably different from a normal 
(Gaussian) distribution because of a tail at positive Sz mean . This is an edge effect: near the 
edges of the mean travel-time perturbation maps, <5r mea n is predominantly positive and large. 
Indeed, narrow white bands are clearly visible on these maps. These bands are produced by 
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Table 3 Linear Pearson correlation coefficients between the different travel-time maps in the quiet Sun, for 
truncated maps: instead of using the full 256 x 256 pixels, we restrict the calculation to an inner square of 
215 x 215 pixels, thus avoiding the map edges. 


Map 

KD97/GB02 

KD97/GB04 

GB02/GB04 

^ r mean A = 6.2 Mm 

0.550 

0.556 

0.958 

S r mea n A = 30.55 Mm 

0.673 

0.639 

0.964 

<5 T'diff A = 6.2 Mm 

0.887 

0.852 

0.953 

<5t(jiff A = 30.55 Mm 

0.873 

0.835 

0.952 


the algorithm we use to compute the cross-covariances: the center-to-annulus averaging is 
performed in the Fourier domain by applying the convolution theorem. However, because 
the Doppler- velocity data cube is not periodic in the horizontal direction, this creates a spu- 
rious behavior at the edges of the map. Ignoring these map edges produces distributions 
that are very close to Gaussian ones. It also somewhat reduces the correlation coefficients 
KD97/GB02 and KD97/GB04 for the mean travel times, but not for the difference travel 
times (see Table 3). The correlation coefficients for the difference travel times remain high. 
At A = 6.2 Mm, this is probably because the difference maps, unlike the mean travel-time 
maps, include a large-scale regular signal (supergranular flows) and not just noise (caused 
by convection and random realization of solar oscillations). At A = 30.55 Mm, there is no 
real signal, and what we see on the maps is realization noise (or, more exactly, filtered re- 
alization noise, since the data cube was phase-speed filtered). Therefore at this distance we 
are testing how the different methods are measuring the noise. 

We also tested the impact of the location and width of the temporal window on the mean 
and difference travel-time values. Figures 4, 5, and 6 show scatter plots for KD97, GB02, 
and GB04. Figure 7 shows the standard deviation of the mean and difference travel-time 
perturbations as a function of the window center and width. 

It appears that at short distance (A = 6.2 Mm) the travel-time definitions are very sensi- 
tive to the temporal window used, especially for the KD97 definition. The standard devia- 
tion strongly depends on this temporal window: we can almost double the noise level on the 
travel-time maps by simply changing the width of the window. At 6.2 Mm, the first-bounce 
skip on the cross-covariances is entangled at short times | / ] with an artifact produced by the 
phase-speed filter (see e.g. Couvidat and Birch, 2009), and it is entangled at larger times 1 t ] 
with the second-bounce skip. The artifact produced by the filter is a group of peaks at small 
t, with the overall shape of a Gabor wavelet and that partly overlap with the first-bounce skip 
peaks, thus affecting their shape. Therefore by shifting the center of the temporal window 
or by changing its width, we more or less select the parts of the first-bounce skip contam- 
inated by this filter artifact (and to a lesser extent by the second-bounce skip). A way of 
alleviating this problem is by using broader phase-speed filters: the filter artifact amplitude 
decreases with wider biters. For a given supergranule, the value of the difference travel-time 
perturbation obtained with GB02 and GB04 can be larger (in absolute value) or smaller than 
the value obtained with the Gabor-wavelet bt depending on how we center the temporal 
window. 

The issue with the use of the phase-speed biter and its interaction with the temporal win- 
dow may potentially be discarded with the HMI data. The phase-speed biter was introduced 
for analysis of MDI data because of the limited spatial resolution of the short distances on 
the cross-covariance diagram. It is expected that the HMI optics are signibcantly improved. 
Thus we may not need to do any phase-speed bitering at short distances (for an example 
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A=6.2 Mm A=30.55 Mm 






Figure 4 Scatter plots for the mean (upper panels) and difference (lower panels) travel-time perturbations 
obtained with the KD97 and GB02 travel-time definitions at two distances A = 6.2 and A = 30.55 Mm. 
The thick solid lines are fifth-order polynomials fitted to the data. The different colors are for different cen- 
ters [ip] and widths [5?] of the temporal window used to select the first-bounce skip during the travel-time 
computation. 


of time - distance diagram obtained with no phase- speed filters, see the one derived from 
Hinode data by Sekii et al, 2007). At A = 30.55 Mm, the filter artifact, first-bounce and 
second-bounce skips, are clearly separated, explaining why the temporal window has less of 
an impact on the travel times. 

The relationships between the KD97 travel times and GB02 or GB04 are not linear. The 
relationships between GB02 and GB04 are linear for small travel times (< five seconds in 
absolute value), as expected since GB04 is a linearized version of GB02. 

4.2. Active Region 

Although GB02 and GB04 were designed for quiet-Sun data, it is possible to use them for 
active-region data, provided that we take into account the decrease in amplitude of the cross- 
covariances that occurs inside sunspots (due to a reduction of the acoustic power). A simple 
way of taking into account this decrease is to normalize the cross-covariances C( A , r) - for 
a given A - at each location r by dividing C(A, r) by its maximum amplitude (L. Gizon, 
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A=6.2 Mm A=30.55 Mm 






Figure 5 Scatter plots for the mean (upper panels) and difference (lower panels) travel-time perturbations 
obtained with the KD97 and GB04 travel-time definitions at two distances A = 6.2 and A = 30.55 Mm. 
The thick solid lines are fifth-order polynomials fitted to the data. The different colors are for different cen- 
ters [tg] and widths [<5r] of the temporal window used to select the first-bounce skip during the travel-time 
computation. 


private communication, 2009): Co(A,r) = C(A, r)/max(C(A, r)). GB02 and GB04 are 
then applied to Co (A , r) instead of C(A , r). The Gabor- wavelet definition is not affected by 
this normalization, because it fits for the amplitude [A] of the wavelet at the same time as the 
travel times. Note that the combination of the decrease in the amplitude of oscillations and 
the use of phase-speed filters may cause systematic shifts of the travel times, particularly, 
the mean travel times (Rajaguru et al, 2006; Parchevsky, Zhao, and Kosovichev, 2008; 
Nigam and Kosovichev, 2010). 

For the active-region test, we use MDI Hi-Res data centered on active region NOAA 
8243 of 18 June 1998. As for the quiet-Sun data cube, the line-of-sight Doppler- velocity 
data cube has the following dimensions: 256 x 256 x 512, with a spatial resolution 5x = 
Sy = 0.826 Mm, and a temporal resolution tit = I minute. The observation duration is 
512 minutes (or 8.5 hours). The reference cross-covariances [C ref ] are obtained by spa- 
tially averaging over the ranges x = [20, 235] and y = [20, 80] in pixel units. These ranges 
are selected because the corresponding MDI map of the line-of-sight magnetic held shows 
no significant magnetic activity there. The maximum value of the cross-covariances at each 
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Figure 6 Scatter plots for the mean (upper panels) and difference (lower panels) travel-time perturbations 
obtained with the GB02 and GB04 travel-time definitions at two distances A = 6.2 and A = 30.55 Mm. 
The thick solid lines are fifth-order polynomials fitted to the data. The different colors are for different cen- 
ters [ip] and widths [5?] of the temporal window used to select the first-bounce skip during the travel-time 
computation. 


location [r] is obtained by Fourier interpolating these cross-covariances on a fine temporal 
grid and then reading their maximum amplitude. Figures 8 and 9 show the mean and dif- 
ference travel-time perturbations obtained with the three definitions. The travel-time maps 
are rather similar especially between KD97 and GB02, showing first a positive mean travel- 
time perturbation at short distances A and then a negative perturbation. Even though KD97 
and GB02 agree well with each other overall, there can be some significant discrepancies 
locally, in some parts of the sunspot. These discrepancies may be caused by the realization 
noise of solar oscillations, which may be higher in sunspots due to the smaller number of 
excitation events in strong field regions, and, perhaps, by some systematic effects since these 
travel-time definitions are not identical. The amplitude of the travel times in the sunspot is 
underestimated by GB04 (in absolute value) compared to the other two definitions. It ap- 
pears that for sunspots, GB02 and GB04 return rather different travel-time perturbations, 
which can be explained by the fact that GB04 strictly reduces to GB02 only when the cross- 
covariances are noiseless and do not have variations in amplitude (GB04 is a linearization 
of GB02 with respect to e for a small e), which is not the case in active regions. In addi- 
tion to the strong non-linear amplitude variations of the cross-covariance in such regions, 
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Figure 7 Standard deviations of the mean and difference travel-time perturbations at A = 6.2 Mm (upper 
panel), and A = 30.55 Mm (lower panel) for five different centers [to] and widths [5f ] of the temporal window 
used to select the first-bounce ridge. 


the noise level of C(A, r) is actually significantly higher than in the quiet Sun. Couvidat 
and Birch (2009) showed that for strong horizontal- flow velocities, the relationship between 
flow strength and resulting travel-time perturbation is not linear, and this non-linearity sets 
in much earlier (for slower flows) for GB04 than for GB02 and KD97. Therefore, we ul- 
timately decided not to use GB04 in the pipeline. Of course, this does not preclude users 
of HMI data from using this or any other procedure in their work. For the linear inversions 
of the mean travel times in the HMI pipeline, we will use both the ray-path (Kosovichev 
and Duvall, 1997) and Born-approximation kernels (Birch, Kosovichev, and Duvall, 2004). 
Since the data-analysis pipeline will also provide the travel-time measurements described in 
this paper, the HMI data users will be able to develop their own inversion techniques. 

Figure 10 shows the standard deviations of the mean and difference travel-time pertur- 
bations for the three travel-time definitions and for the 1 1 distances [A] studied here. These 
standard deviations are computed in a quiet-Sun region surrounding, but excluding, the ac- 
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Figure 8 Difference travel-time perturbations for the active region NOAA 8243 obtained with KD97 (left 
column), GB02 (central column), and GB04 (right panel), and for five source - receiver distances [A]. The 
colorscale is the same for all the panels, and ranges from —40 seconds to +40 seconds. 


tive region. They confirm the conclusion of the previous section in which we showed that 
at short distances [A] the travel-time maps produced by the Gabor-wavelet definition have 
a smaller standard deviation than those produced by GB02 and GB04, while the opposite is 
true at larger A. 
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Figure 9 Mean travel-time perturbations for the active region NOAA 8243 obtained with KD97 (left col- 
umn), GB02 (central column), and GB04 (right panel), and for five source -receiver distances [A]. The 
colorscale is the same for all the panels, and ranges from —40 seconds to +40 seconds. 


5. Conclusion 

We have implemented three travel-time measurement procedures (Kosovichev and Du- 
vall, 1997; Gizon and Birch 2002, 2004). The subroutines computing these travel times of 
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Figure 10 Upper panel: standard deviation of the mean travel-time perturbations for KD97 (solid line), 
GB02 (dashed line), and GB04 (dash-dotted line). Lower panel: same but for the difference travel-time per- 
turbations. 


wavepackets are part of the HMI time - distance helioseismology pipeline at the SDO Joint 
Science Operations Center at Stanford University. Using quiet-Sun and active -region data 
from the SOHO/MDI instmment we compared the mean and difference travel-time pertur- 
bations obtained by using these three definitions. They all generate similar travel times in 
the quiet Sun, provided that, at short distances [A], we use an appropriate temporal window 
to isolate the first-bounce skip on the cross-covariances. Indeed, at A = 6.2 Mm, the travel 
times strongly depend on the location and width of this temporal window. This is mainly due 
to the phase-speed filter applied to the data cube. This phase-speed filter is used to isolate the 
wavefront signal from the artifact signals of the MDI instrument (Duvall et ai, 1997). With- 
out such filter, the travel-time measurements at short distances are not possible from MDI 
data. However, this filter introduces issues of its own: the travel times measured at short dis- 
tances depend on the parameters of these phase-speed filters. Using broader filters than those 
considered here reduces the magnitude of this problem. HMI data might not necessitate the 
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application of phase-speed filtering, which should give us access to more robust travel-time 
estimates at short distances. The distributions of the mean and difference travel times ob- 
tained in the quiet Sun are close to a Gaussian distribution (especially when the map edges 
are discarded) and are very similar for the three travel-time definitions. The Gabor-wavelet 
technique (Kosovichev and Duvall, 1997) produces travel-time maps with a standard devia- 
tion smaller than the Gizon and Birch (2002, 2004) definitions for A < 15 Mm, but with a 
larger standard deviation for larger A. The relationships between KD97 and GB02 or GB04 
are not linear, while the relationships between GB02 and GB04 are linear for (absolute) 
times smaller than about five seconds. In the active region studied here, GB04 returned 
travel times that are significantly smaller (in absolute value) than those of GB02 and KD97. 
We ultimately decided not to include GB04 in the pipeline. 
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